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Abstract 

The effects of dissipation on tlie space time evolution of matter formed in ultra- 
relativistic lieavy ion collision is discussed. The thermal photon spectra for RHIC 
and LHC energies with viscous flow is considered. The effects of viscosity in the 
thermal single photon spectra is seen to be important in QGP phase as compared to 
the hadronic phase. Recently available data from WA80 group at SPS energies for 
S + Au collision are compared with theoretical calculations. The experimental data 
do not appear to be compatible with the formation of matter in the pure hadronic 
phase. 
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I. Introduction 



The primary motivation for studying heavy ion colhsions at ultrarelativistic energies 
is that such reactions provide an unique opportunity to probe hadronic matter at 
very high temperature and/or density. Numerical simulations of the QCD equation 
of state on the lattice predict that at a sufficiently high temperature (~ 160 — 200) 
MeV, a phase transition should occur from the colour confined chirally asymmetric 
phase of QCD, the hadronic phase, to a (locally) colour deconfined, chirally sym- 
metric phase, referred to as the Quark-Gluon Plasma(QGP) The QCD phase 
transition is of great relevance in cosmology also [Q; within the standard big bang 
model, the early universe should have experienced such a transition a few microsec- 
onds after the big bang. Although the order of this phase transition remains to be 
settled, a substantial volume of work has already been carried out, assuming it to 
be of first order both in the laboratory as well as in the early universe. In this work, 
we shall tacitly assume that the hadron to QGP phase transition is indeed of first 
order. 

In most of the works aimed at devising signals of the above mentioned phase 
transition in the laboratory, the usual scenario is as follows: after a canonical 
(proper)time ~ 1 fm/c (the strong interaction time scale), the excited matter within 
the reaction volume of the colliding heavy ions comes to a (local) thermal equilibrium 
with a temperature Tj. If this temperature Tj is larger than the critical temperature 
Tc of the QCD phase transition, then the excited matter would exist in the form of 
QGP. The QGP then evolves in space and time, cooling down to the critical tem- 
perature Tc when the confinement phase transition starts; for a first order transition 
the system remains in a mixture of QGP and hadronic matter at a temperature 
for some time the duration of which is determined by the relative number of the 
degrees of freedom in the QGP and hadronic phases. Such a configuration continues 
until all the latent heat which maintains the temperature at by compensating the 
cooling due to expansion is fully exhausted. At this stage all of the QGP converts 
to hadronic phase and the system continues to cool further, its dynamics being gov- 
erned by the hadronic equation of state, until the mean free paths of the constituent 
hadrons become too long to maintain a collective behaviour. At this point, usually 
referred to as freeze-out era, the momentum distributions of the particles become 
frozen at the characteristic values when the particles suffer the last collision, and the 
particles free stream thereafter, carrying the space time integrated informations to 
the detectors. For the sake of simplicity, most studies on QGP diagnostics assume 
that the system behaves as an ideal i.e. nondissipative, fluid throughout the entire 
evolution. 

Given the tremendous importance of the issue, it is imperative to understand 
how far, if at all, such an idealisation is justifled. Even though a rigorous treatment 
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of the dissipative effects in relativistic fluid dynamics is still beset with several 
technical as well as conceptual difficulties ^, one must admit that in any realistic 
scenario, dissipative effects, in principle should have an important role. In this work 
we consider the relativistic fluid to be weakly dissipative, as has been argued to be 
the case by earlier authors]^ and explore the ensuing consequences on the estimation 
of the initial temperature and the photon spectra, data for which have become 
available of late. Since our motivation in the present work is to study the relative 
importance of the dissipative effect vis-a-vis an ideal fluid dynamic scenario that 
has been routinely used by a number of authors 0, |^ |], |^ to compare the currently 
available data and/or to estimate experimentally measurable quantities at RHIC or 
LHC, we adopt the same evolution scenario (Bjorken scaling flow), together with 
the same formation time Tj (~ 1 fm/c). It must be mentioned at the outset that 
our present purpose is demonstrative rather than making firm phenomenological 
comparisons. 

The organisation of the present work is as follows. In the next section, we 
briefly recapitulate the space-time evolution of the system with dissipation. In 
section III we discuss the effect of dissipation on the boundary or initial conditions 
governing the evolution equation. Section IV briefly recapitulates the formalism for 
photon diagnostics of QGP. Finally in section V, we present the results of our actual 
estimates at LHC, RHIC and SPS energies and compare them with the available 



datafT^ at SPS. Section VI contains a brief summary and conclusions. 



II. Viscous Hydrodynamics at Relativistic Energies 

For a perfect fluid, the mean free path of the constituent particles is negligibly small 
compared to all the available length scales in the system. In a finite system like 
colliding nuclei, however large, such a criterion may not be fulfilled and as such, the 
changes in the collective quantities like pressure, energy density, number density, 
velocity etc. over a distance of a mean free path cannot always be ignored. This 
should lead to an inherent dissipative behaviour. 

In the presence of dissipation, the kinetic energy of the fluid decays as heat 
energy. This requires a redefinition of the energy-momentum tensor and the particle 
current [|l^. For the sake of brevity, we simply quote the result in the form of 
energy-momentum conservation in an imperfect fluid which reads 

Eq.(l) is the scaling Navier-Stokes(NS) equation where we have assumed the validity 



of the scaling solution U^, usually referred to as the Bjorken solution |13|. The rate 



of heat generation is (4r7/3 + C)/t^, which vanishes in the case of an ideal fluid 
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(?7 ^ 0,C — ^ 0). Obviously, the occurrence of the viscosity co-efficients rj and C in 
eq.(l) violates the invariance under r — > — r, and consequently the reversibility [|I^ . 
As a result, entropy is generated during the temporal evolution of the system; thus, 
estimating the initial temperature from the final multiplicity in the usual manner 
[|] constitutes an overestimate (see Appendix-I). 

We assume, without any loss of generality that the functional relation between 
the pressure and energy density in case of dissipative flow is same as that of ideal 



flow [11|. We shall, as is the usual practice, assume that the baryon density is zero 



in the central region so that the thermal conductivity i^' = throughout. 

Ila. Quark-Gluon Plasma Phase 

The limits on the shear viscosity for the NS equation (eq. 1) in the QGP phase have 
been discussed in depth in [§](see also JT^). The acceptable range of r/ is 



2T^ <T]< 3T%tT) (2) 

Obviously, the range is quite small unless the product Tr (very late times and/or 
large temperatures) is very large. As should be clear from what follows, such a 
situation is rather unlikely within the QGP phase. 

The bulk viscosity (, in general can be as large as rj, or in certain circumstances, 
even larger [|1^. It was however shown by Weinberg some time ago that if the 



trace of the energy momentum tensor be expressible as a function of e and the 
pressure is the same function of e as in the adiabatic case, then C = 0. We thus 
obtain for the scaling NS equation in the QGP phase 

de e + P 4 2 /q\ 

Although eq. (3) is somewhat simpler than eq. (1), we work with the more general 
form for the time being. For e = SgqT^ and P = e/3, where ag = (vr^/90)(7QGP, we 
can solve eq. (1) with the boundary condition T = Ti sX t = Ti. 



T = Un/rf' + i^VQo + kgo)-- - - (4) 

D O an Ti 



r J 

where tjqq = rig/T^ and (qq = (q/T^. The proper time tq taken by the system to 
cool down from the initial temperature Ti{> T^) to the critical temperature Tc can 
be calculated by inverting eq. (4): 

'Tq = 'rQ{Ti,T,,T^) (5) 

The life time of the QGP phase is tq—Ti. It can be readily seen that in the presence of 
dissipation, the life time of the QGP phase becomes longer (due to viscous heating) 
relative to the ideal case, for the same values of Tj, Ti and Tc. For an actual estimate 
we use rjQo = 2.5 , as dictated by eq. (2) and Cqo = in the pure QGP phase. 



4 



lib. The Mixed Phase 



When the temperature of the system comes down to T^, the system goes over to a 
mixed phase composed of QGP and hadronic matter. In this phase the temperature 
remains constant but the energy density changes due to the change in the statistical 
degeneracy on account of hadronization. We define fq as the volume fraction of the 
QGP in the mixed phase, so that the energy density at any proper time r during 
the life time of the mixed phase is given by 

e(r) = /Q(r)eQ + (1 - /Q(r))eH (6) 

where eQ^en) is the energy density in the QGP(hadronic) phase at temperature 
Tc. Eq.(6) is obtained by a linear interpolation between the two points (l,eQ) and 
(0,6/^) in the (/Q,e) plane. For the coefficients of viscosity, we can formally use a 
similar interpolation, 

Vir) = fQiT)r^Q + (l - fQ{T))rj„ (7) 

and 

C(r) = /Q(r)CQ + (1 - /Q(r))CH (8) 

where the subscript Q{H) denotes the value in the QGP(Hadronic) phase at Tc. 
The applicable limits on tiqq and Cqo have already been discussed in the previous 
section. The corresponding situation in the hadronic phase is explained in the 
following section. 

The evolution equation for /q(t) in the mixed phase can be readily obtained 
from eqs. (1), (6), (7) and (8): 

where a = (en + Pe)/Ae ; b = A^riQ - r/j^)/3Ae ; c = (2^0 + |r?//)/Ae ; Ae = eg - en- 
The solution of this equation with the boundary condition /q = 1 at r = tq is given 
by 

/q = {ab - c) \Ei{b/T) - Eiib/TQ)] ^ + (a + lf-9.ei^/-Q~^/-) - a (10) 

r r 

where Ei{x) is the exponential function. For 6 = c = 0(?7 = ^ = 0) we recover the 
solution for fq in case of a perfect fluid. Note however that the bulk viscosity 
( in the mixed phase may be non-zero due to the change in sound velocity over a 
finite relaxation time |T^. We use ( = 2ri/3, following Hosoya and Kajantie [|15 



The time th, when the mixed phase ends, is obtained from eq. (10) by demanding 
fQ^Tn) = 0. We thus have, 

TH = rniTq) (11) 

for given equations of state in the QGP and the hadronic phases. The life time of 
the mixed phase is then th — tq. 
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lie. Hadronic Phase 



The pure hadronic phase starts at T = at the time tjj so that the solution of eq. 
(1) in the hadronic phase is again eq. (4) with the appropriate substitution of the 
ag by an, Ti by and n by th: 



1 1 . . 1 1 



T = T,(rH/r)i/3 + (^r/H0 + K. 



D 8 clh Th 



T T 



(12) 



In the hadronic phase, the co-efficient of the shear viscosity can be obtained from 
the estimated transport cross-section according to the relation rjn ^ T/atr ~ 
(T/200Me\/) (0.5 - l)/fm3||]. The expected range of atr ~ 10-20 mb which 
corresponds to rjuo ~ 0.75 — 1.5 at Tc. In what follows we shall use the larger value 
of 1.5 for 7] HQ as a conservative case. The value of the bulk viscosity can be non-zero 
in the hadronic phase due to lack of chemical equlibrium and its consequence on the 
sound velocity [|T8], we again take C,h = '^rjn/'i iT5[ . 



As mentioned in Section I, the system freezes out at a time denoted by r/, when 
it becomes too dilute to maintain a collective flow. There is as yet no universal 
criterion for determining the instant of freeze-out, for a detailed discussion, see [^. 
In the present context, we describe freeze-out through a freeze-out temperature Tf 
which we treat as a free parameter (< 140 MeV). Then from eq.(12) Tf is given by 

= yi^H,Tc,Tf) (13) 

so that the life time of the hadronic phase is Tf — th- 



III. Hydrodynamic Flow and Connections with Ob- 
servables 

In the case of an ideal fluid, the conservation of entropy implies that the rapidity 



density dN/dy is a constant of motion for the isentropic flow [|T3|. In such cir- 
cumstances, the experimentally observed multiplicity, dN/dy may be related to a 
combination of the initial temperature T, and the initial time as T^Ti. Assuming 
an appropriate value of (taken to be ~ 1 fm/c), one can estimate Tj. 

For dissipative systems, such an estimate is obviously inapplicable (see appendix 
I). Generation of entropy during the evolution invalidates the role of dN/dy as a 
handy constant of motion. Moreover, the irreversibility arising out of dissipative 
effects implies that estimation of the initial temperature from the flnal rapidity 
density is no longer a trivial task. We can, nevertheless, relate the experimental 
dN/ dy to the freeze-out temperature Tf and the freeze-out time Tf by the relation, 

^ = TTRlAanTfTf/c (14) 
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where Ra is the radius of the colhding nuclei (we consider A A colhsion for simphcity) 
and c is a constant ~ 3.6. The parameter qh occurring in eq. (14) is the so called 
statistical degeneracy (a = (7r^/90)(?) for an effectively massless hadronic gas. In 
this work we assume the hadronic matter to consist of vr, p, uj and r/, which allows 
us to estimate an through a parametrization introduced earlier [|19|. It must also 
be mentioned that eq. (14) still constitutes a tremendous amount of idealization; 
for a discussion of all the caveats, see0. For our present purpose, we assume that 
eq. (14) does hold at the freeze-out era. We shall also neglect transverse expansion 
in this work, which would amount to an overestimation of rj, for a given Tj and 
dN/dy. We however do not consider this to be a major factor; this is to be discussed 
in some detail in Sec. V. 



IV. Photon Spectra 



In this section we will briefly discuss the photon spectra due to an expanding quark 
gluon plasma. In the QGP phase the main contribution comes from the annihila- 
tion {qq — * (77) and Compton processes {q{q)g In the hadronic phase 
(composed of vr, p, u, 77), an array of reactions like vrp — > TC'y, tti] — > 717, tttt 777, 
TTTT — )■ p7, TTTT —>■ 77 and the decays p° — > 7r"^7r~7 and uj ir^'-f are considered pO| . 
The rate of emission of photons per unit volume from QGP phase is given by 




d^p 9 27r2^ ^ g' ^ ' ^ ' ^^^^ 

For our exploratory calculations we will consider the rates of photon emission from 
the QGP and the hadronic matter to be equal as in pO, It is a very useful first 



approximation for identifying the pt window where we concentrate on the signature 
of QGP. The pt distribution of photons is obtained by convoluting the basic rates 
with the space time history. 



dN 



d^PTdy 



TTRi 





/o) e(M) 



rdr drj 



(16) 



where 9(M) = 0(eQ — e)0(e — en)- The above equation along with eqs. (4), (10) 
and (12) is used to evaluate the photon spectra. 
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V. Results 



As we have mentioned in section III, estimation of the initial energy density /tem- 
perature from the final state rapidity density is no longer a trivial task in the 
presence of dissipation. The principle however is quite straight forward. In our 
algorithm, we tacitly assume that the initial formation time of the thermalised 
system , be it in the QGP phase or the hadronic phase, has the canonical value of 
1 fm/c. That the initial formation time (and consequently the initial temperature 
Tj) is a poorly known quantity is beyond dispute. Indeed, it is now widely believed 
[ p^ , p3| , that at RHIC or LHC, the gluons may equilibrate considerably earher 
than 1 fm/c. However, even at these energies, the u/d quarks appear to equilibrate 
only ~ 1 fm/c [^]. Moreover, all the works aimed at comparing the experimental 
data with theoretical predictions take, to be, as a rule, equal to or greater than 1 



fm |g, ^, 1^, |T3|, |2^. Most of these works also use the Bjorken scaling hydrodynamics 
to emulate the space-time evolution of the system. For a meaningful comparison 
with these works, we should also therefore adopt the same scenario Q. 

We treat Tj as a parameter; for each Tj, we let the system evolve forward in 
time under the condition of dissipative fluid dynamics (eq. (1)) till a given freeze- 
out temperature Ty is reached. Thus tj is determined. We then compute dN/dy 
at this instant of time from eq. (14) and compare it with the experimental dN/dy. 
The value of Tj for which the calculated dN/ dy matches the experimental number 
is taken to be the approximate initial temperature. Schematically the calculation 
for Tj proceeds (forward in time) as follows. 

TQiT,) ^ MTQiT,)) - TfiTHiTQm) (17) 

i.e. th depends on Tj through tq and finally Tf depends on Tj through th {cf eqs. 
(5), (11) and (13)). We then solve the non-linear eq. (14) to estimate Tj for given 
values of dN/dy and Tf. We adopt a truely conservative view point in the sense 
that we start with values of Tj less than T^ so that the system is formed in the 
hadronic phase, and allow values of Tj > Tc only if the discrepancy between the 
calculated and experimental (or estimated) values of dN/dy necessarily demands it. 
In these circumtances, we naturaly let the system evolve first in the quark gluon 
plasma phase, then the mixed phase and the hadronic phase. 

The result of these calculations are summarized in table I. For the sake of com- 
parison, we also list the corresponding values as calculated in the ideal case (su- 
perscript Bj). For RHIC and LHC , the experimental dN/dy values are the usual 

^ The assumption of 1 fm/c as the initial formation time need not be taken seriously. For our 
present purpose, one can always adopt the viewpoint that we compare the temperature at a proper 
time 1 fm/c, a sufficiently early instant of time in the evolution history, in the dissipative vis-a-vis 
ideal fluid scenarios. 
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estimates obtained from extrapolating the pp data||2^. For SPS , the dN/dy is 
known from actual measurement [P7|] . 

It can be readily seen from table I that in all cases, the initial temperatures for 
the dissipative case are lower than in the ideal case. Lifetimes for the dissipative 
case are also comparatively longer. One also finds that the lower the freeze-out 
temperature, the lower the initial temperature although the dependence is rather 
weak. 

It is evident from Table I that both in RHIC and LHC, the initial temperature 
is substantially higher than even for dissipative dynamics. (We use the canonical 
value of 160 MeV for in this work.) 

The situation at SPS demands special discussion. Please therefore refer to Table 
2a. We find that an initial QGP phase may be formed only if we use Tc lower than 160 
MeV. For example, T^. of 150 MeV allows a QGP phase with an initial temperature 
of 154 MeV only, and that too for a freeze-out temperature of 100 MeV. Thus an 
unambiguous conclusion about the formation of the QGP on the basis of presently 
available SPS data appears to be doubtful. We therefore explore the situation in 
somewhat greater detail in Table 2b. 

In Table 2b we study, on the one hand, the possibility that the system is formed 
entirely in the hadronic phase with dissipative dynamics. We look at the most con- 
servative scenario of very low freeze-out temperature, consistent with the argument 
of Hama and NavarrapSj for the case of S* + Au system. We find that to obtain 
a dN/dy of 225, an initial temperature of 262 MeV is required which is obviously 
much larger than the usual Tc of 160 MeV (or even 200 MeV ). It is thus fair to 
conclude that for a reasonable formation time of 1 fm/c , the present SPS data are 
not compatible with a purely hadronic phase even if dissipative effects are taken into 
account. This, together with the results detailed in Table 2a makes the situation 
most interesting indeed. We thus explore the remaining possibility that the initial 
system is formed in the mixed phase, i. e., Ti = T^. 

Table 2b shows that this case fits the data reasonably well. Depending on the 
choice of Tj, the initial system could be a mixture of 80 — 90% quark matter and 
the rest hadronic matter, co-existing in thermal equilibrium at a temperature T^. It 
thus appears that there is a strong hint of deconfinement phase transition of first 
order in the currently available data but for a clear indication of the formation of the 
QGP phase , one must go to higher energies at RHIC or LHC. This situation is also 
similar to that advocated by Werner [p9|] of late, where he suggests the formation of 
quark blobs in hadronic matter at SPS energies. 

Effects of these considerations on QGP diagnostics are of direct relevance in the 
present scenario. One of the promising signals of QGP formation is direct photons 
0. We thus show in figs. 1 and 2 the calculated photon spectra for ideal as 
well as dissipative fluid dynamics at RHIC and LHC energies respectively. The 
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immediate feature that stands out from these figures is that the effect of dissipation 
is stronger in the QGP phase than in the hadronic phase. This is to be expected as 
the temperature in the QGP phase is higher, corresponding to the larger values of 
the co-efficients of viscosity. It should also be noted that even if the estimate of Tj 
is lower than the ideal case by only 10% or less, the change in the photon spectra 
at high pt is quite sizeable. 

In fig. 3, we compare the photon spectra corresponding to the configuration of 
Table 2b (initial mixed phase at a temperature of 160 MeV ) with the experimental 
data of WA80|T0|, 0. For the sake of comparison, we also show the photon spectra 
from purely hadronic phase with the initial configuration of Table 2b. Obviously, 
the agreement with the initial mixed phase is better, considering the fact that the 
experimental data gives the upper limit for the photon spectra (fig. 3). Before closing 
this section, we feel that afew comments on our neglect of transverse expansion are 
in order. As we already mentioned, this amounts to an overestimate of rj, especially 
for RHIC and LHC energies. It is well known that the transverse flow effects at SPS 
energies are not at all large |21|. Nevertheless, an inspection of Table 1. reveals 
that the change due to dissipative effects in total lifetime or Tj relative to the ideal 
case is less than 10% at both RHIC and LHC energies. Since the same effect is seen 
both in Tf and th-, it is reasonable to expect that even in the presence of transverse 
flow, the same qualitative features would obtain. This however ought to be verified 
in a full (3+1) dimensional calculation which we plan to perform in the near future. 



VI. Summary 

In this work, we have explored the effects of dissipation in the space time evolution 
of matter formed in very energetic collisions of heavy ions. We find that even with 
the assumption of scaling hydrodynamic flow and weak dissipation, the deviation 
from the ideal case is not negligible. The presence of dissipation reduces the rate of 
cooling, resulting in a longer total lifetime. For a given total multiplicity in the final 
state, the estimated initial temperature for the dissipative case is smaller than in 
the ideal situation. This effect is most pronounced at SPS energies. It appears that 
the presence of dissipation at RHIC or LHC energies may not change the estimated 
Ti to such an extent as to qualitatively invalidate the conclusions derived from the 
assumption of ideal hydrodynamics. At the level of the present analysis, we find that 
the measured dN/dy of 225 cannot be easily reproduced without a first order phase 
transition (Table 2b). The most likely situation seems to be an initial configuration 
of a mixed state of hadronic and quark matter. 

In the tharmal single photon spectra, we observe that the effect of dissipation 
should show up at the high pt sector. This is quite natural, given the fact that 
dissipation results in a lower initial temperature. In all cases, the effect of viscosity is 
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important in the QGP phase as compared to the hadronic phase. Recently available 
data from the WA80 group at SPS energies for S + Au collisions have been compared 
with theoretical calculations without any additional free parameters, we find that 
the parameters extracted from the dN/dy analysis (Table 2b) yield a reasonably 
good fit to the data, if the initial configuration is indeed a mixed phase. The 
data are not compatible with formation of matter in a pure hadronic phase. These 
conclusions arc not expected to be materially alterred due to transverse fiow which 
we have not taken into account for the present. 
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Table Captions 



• Table 1: Comparison of time scales for the different phases and initial temper- 
atures between Bjorken hydrodynamics(superscript Bj) and viscous hydrody- 
namics at RHIC and LHC energies. dN/dy is the multiplicity, Tj, Tc and Tf 
are the initial, critical and freeze-out temperatures respectively. rj,rQ, th and 
Tf are the initial thermalisation time, time when phase transition starts, time 
when mixed phase ends and freeze-out time rescpectively. 

• Table 2a: Same as Table 1 at SPS energies. The value of dN/dy is taken from 
experiment. 

• Table 2b: Same as Table 2a when matter is formed in pure hadronic or mixed 
phase. 



Figure Captions 

• Figure 1: Photon spectrum at RHIC energies for ideal and viscous 
flow. The inputs are taken from table. 1, with freeze-out temperature Tf = 
100 MeV. QM(HM) denotes the contributions from pure QGP(hadronic) 
and QGP (hadronic) part of the mixed phase . 

• Figure 2: Photon spectra at LHC energies with and without viscous flow. 



Figure 3: Photon Spectra at SPS energies with viscous flow and comparison 
with WA80 data. The dashed curve denotes the photon spectrum with inputs 
from no phase transition scenarion of table 2b. The solid curve denotes the 
photon spectrum with inputs from table 2b. The freeze-put temperature is 
taken as 65 MeV as calculated by applying the method of ref . . The arrow 
denotes the experimental data from WA80pO|, p!0| ]. 
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APPENDIX-I 



In case of ideal fluid, the conservation of entropy implies that the rapidity density 
dN/ dy is a constant of motion. In such cases the experimentally measured value of 
the final rapidity density, which is a quantity depends on the integration over the 
space time history, can be connected to the initial temperature. In these calculations 
one assumes the validity of the extrapolation from the freeze point (with temperature 
Tf and time Tf) to the initial point spccficd by temperature Tj and proper time 
backward in time. But in the presence of dissipative effects the process is irreversible, 
and one shoud take care of this fact. In this appendix we will show that if we solve 
the eq.(l) backward in time then we always overestimate the initial temperature for 
a given vaule of dN/dy and Tj. 

Under the transformation r — r, eq.(l) becomes, 

de (e + P) ,4 ^, , n , , 

For the equation of state of a massless gas, eq.(15) can be written as 

where 77 = ^0^^, C = Co^^, b = |?7o + Co, ^fc = (7r^/90)(?fc, gk is the statistical 
degeneracy. The solution of eq.(A.2) with the boundary condition T — Tf at t — Tf 
can be written as, 

Tr'/' = K {dN/dy - — (rf^' - r-'/') {A.3) 

where we have used the relation TfT^J^ — K{dN j dyY^^ ■ The initial temperature 
estimated from eq.(A.3) is 

Tlrl'^ = K {dN/dyf^ - ±- {rf' - rf'^) {AA) 

T/ denotes the initial temperature when we solve the hydrodynamic equation back- 
ward in time. Similarly we can solve the hydrodynamic eq.(l) with boundary con- 
dition T = Tj at r = Tj, the solution is given by. 

At the freeze-out time eq.(A.5) can be written in terms of dN/dy, 

T^" = K {dN/dy + ^ {rf' - rr^/') (A.6) 

In case of ideal hydrodynamics, (6 = 0) T/ = Tj. But in case of viscous flow for a 
given dN/dy and Tj 

t; > Ti {A.l) 
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Table 1 





(MeV) 


T 

-'- c 

(MeV) 


(fm/c) 


(fm/c) 


(fm/c) 


(fm/c) 


t./t;^^' 

(MeV) 


RHIC 
(f = 1735) 


100 
140 


160 
160 


1/1 
1/1 


3.5/4.2 
3.6/4.2 


35/34 
36/34 


150/138.4 
54.7/50.5 


225/258 
228/258 


LHC 
(f = 5624) 


100 
140 


160 
160 


1/1 
1/1 


13.7/13.6 
13.8/13.6 


117.2/109.8 
118.2/109.8 


486.8/449.8 
177.4/163.9 


358/382 
359/382 
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Table 2a 



(MeV) 


T 
(MeV) 


(fm/c) 


(fm/c) 


(fm/c) 


(fm/c) 


(MeV) 


SPS 100 
= 225) 100 


120 
150 


1/1 
1/1 


3.5/1.7 
1.1/2.5 


37.9/13.6 
18/20 


67.5/23.5 
67.5/67.5 


165/203 
154/203 



Table 2b 



No 


Phase 


Transition 










dN 
dy 


T 
(MeV) 


n 

(fm/c) 


Tf 
(MeV) 


^f 
(fm/c) 






225 


264 


1 


65 


246 






Formation 


in 


Mixed 


Phase 








dN 
dy 


Tf 
(MeV) 


T,{= T,) 
(MeV) 


(fm/c) 


th 
(fm/c) 


^/ 
(fm/c) 


/o 


225 
225 
225 


140 

100 
65 


160 

160 
160 


1 
1 
1 


16 

15 
14 


25 

68 
247 


0.92 

0.81 
0.75 
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1 2 3 4 5 6 

Pt (GeV) 




Pt (GeV) 



10 



10 



10 



10 



10 



10 



10 



10 



10 







-1 



-2 



-3 



-4 



-5 
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no phase transition 

formed in mixed phase 



WA80 data ( S + Au) , 
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Viscous flow 
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